Collective modes as a probe of the equation of state for partially polarized Fermi gases 
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We calculate the collective modes of a partially polarized Fermi gas trapped in a spherically 
symmetric harmonic potential. We show that the breathing mode frequency exhibits non-monotonic 
dependence on polarization in the entire BCS-BEC crossover region. Moreover, we find that the 
breathing mode can be used to distinguish between two commonly used unitary gas equations of 
state. 
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Recent experimental studies of partially polarized 
Fermi gases P, 01 at ultra-cold temperatures near a 
Feshbach resonance provide a unique window into the 
behavior of strongly interacting fermions. Relying on 
long spin relaxation times, experimentalists polarize their 
two-component Fermi gasses, driving them from super- 
fluid to normal. Consistent with theoretical predictions 
0, 0, 0, 0, Q , they see the cloud phase separate into con- 
centric shells, however, the detailed structure of the shells 
are not completely understood. The Rice experiments 
have largely been analyzed in terms of a two-shell model 
where an unpolarized superfluid core is surrounded by 
a completely polarized normal shell. The MIT exper- 
iments strongly suggest a three-shell structure with 
a partially polarized normal shell between the other re- 
gions. There currently exists no reliable model for the 
equation of state of these partially polarized gases. Here 
we explore to what extent collective modes can be used 
to experimentally probe this equation of state. 

We use a hydrodynamic approach to calculate the 
breathing modes of a partially polarized Fermi gas at 
unitarity. We work with the two most commonly used 
equations of state finding a significant (~ 5%) 

difference between the oscillation frequencies. This sen- 
sitivity should be contrasted to the case of unpolarized 
gases where the differences between the predictions of 
mean field theories and quantum monte-carlo are much 
smaller (~ 2%) Ongoing experimental efforts 

which are attempting to observe these smaller beyond 
mean- field effects 0, could easily distinguish be- 
tween the two equations of state which we consider, and 
can help refine our theories of resonant fermions. 

We study both the polarization and interaction 
strength dependance of the breathing mode frequency. 
We find that the breathing mode frequency is a non- 
monotonic function of polarization for all interaction 
strengths. 

We consider a zero temperature gas of fermionic atoms 
of mass 771 in two hyperfine states \a =t,i> confined 
by a spherically symmetric harmonic potential U{r) = 
[mujQ 1 2)r'^ , where is the trapping frequency. In terms 
of the numbers of atoms N„ in each hyperfine state, the 
polarization is defined as P = {N^ - N^) / {N^ + N^) > Q. 
Assuming local equilibrium, the dynamics of this system 



will be described by a continuity equation dtPcr+'^ ■ja- = 
and an Euler equation mdt^cr + (m/2)Vw^ -I- Vt/eg = 0, 
where pcr,^cr, and jcr = Pa^a are the mass density, veloc- 
ity, and mass current of the atoms in state a. In the local 
density approximation L/cff — U{r) — pcr- Linearizing the 
equations around equilibrium density, pa = p'^ + Spa, and 
writing the density fluctuations in terms of fluctuations 
in chemical potential, Spa- = ^^{dpa / dpu)^Pv, we have 



' Kcr^d^Spi, = V 



— V6pa 

m 



(1) 



where we have introduced the compressibility matrix 
Kst — dps/dfif Explicitly considering the case of har- 
monic trapping, the equilibrium local chemical potentials 
for up and down atoms in the local density approxima- 
tion (LDA) are given by p^{r) = po + h— (ma;§/2)r^ and 
Pl{r) = pq — h — (mu;Q/2)r^. If we know the equation of 
state, we can extract the local equilibrium densities from 
these local chemical potentials. Previous studies have 
shown that this LDA approach is an excellent approx- 
imation unless the number of particles is small and the 
trap deviates significantly from spherical. Under such cir- 
cumstances, surface tension between the superfluid and 
normal regions must be considered 0. Here we neglect 
such finite size effects. The parameters pQ and h are de- 
termined from a constraint on the total number of atoms 
N and the polarization P. 

We formally reduce Eq.Q to a dimensionless form 
by defining, pa- = ■m{mp)^^'^ f„{ p/e, h/ e), = 
iri^lmpy/^ g^uip/^, h/e), and x — \J muj^jlpisr , where 
ja and are dimensionless functions. The average 
chemical potential p = -I- p\)l'i — Mo(l — x^) and 
e — h? /mal with s-wave scattering length and Eq. 
(P) becomes 
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We see that dimensionless oscillation frequencies uj/ljo 
only depend on /^o/e and h/t, or equivalently 



2mpo/h^as and P. 
We work with two different equations of state: a BCS 
mean- field result, and a semi-empirical model introduced 
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by Frederic Chevy [5j . The latter equation of state is suf- 
ficiently simple that we can write the solutions of Eq. Q 
in terms of hypergeometric functions. The BCS mean 
field theory is sufficiently complicated that we must solve 
Eq. ^ numerically. By using spherical symmetry we re- 
duce Eq. iQJ to a differential equation for the radial vari- 
ation of After discretizing space and fourier trans- 
forming with respect to time, this radial equation has 
the form of a matrix eigenvalue problem. Using standard 
sparse matrix techniques we extract the eigenvalues and 
eigenvectors. We verify that our results are independent 
of the discretization procedure. 

Neglectin g th e possibility of a modulated order param- 
eter (FFLO)|]J,the mean-field equation of state is found 
by numerically solving the following equations for the 
density p = pi + pi, the density difference pd = p^ ~ pi, 
and the gap A, 
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In these equations, Ek{r) 
£.ah + ^/{ek-pY 



= {Ek^+Eki)/2, with Eka = 
where tk — h'^k'^/2m, ^| — 



-A2, 

-pi)/2 andh — {pi~pi)/2. The mo- 
menta of the Fermi surfaces are k±{r) = {zLyh^--S^ + 
pY^"^. In the normal state, A = 0, this approach re- 
produces the equation of state of a noninteracting Fermi 
gas. At unitarity, this theory predicts a three-shell struc- 
ture, similar to that observed at MIT however the 
intermediate partially polarized normal shell is predicted 
to be extremely small, and the mean-field theory effec- 
tively drops all interactions in this shell. Figure 1 shows 
the results of this hydrodynamic calculation at unitar- 
ity as — > cx), while figure 2 shows the results at several 
different values of as- 

Our second approximate equation of state was intro- 
duced by Frederic Chevy Unlike the mean-field the- 
ory, this approach assumes that only two phases exist 
at unitarity: a completely unpolarized superfluid, and 
a completely polarized normal state. The forms of the 
equations of state of these two phases are known exactly: 

3/2 

in the normal phase, p = Bp^ j while in the superfluid 
phase, p = A[{p^ + pi)/2Y/'^. Elementary statistical me- 
chanics gives B = (l/67r2)(2rn/?i^)3/2, while dimensional 
analysis requires A — 2Bq^^/'^. Comparison with ex- 
periments and numerical Monte-Carlo calculations give 
« 0.45 Q. In this model, the edge of the superfluid 
corresponds to the edge of the minority cloud (i?x)- 
this boundary the ratio of the chemical potential is a uni- 
versal number, pjp^^ = 2{B/Af/^ - 1 = « -0.061 
0,0. The edge of the majority species cloud {R-\) occurs 



when /i| = 0, allowing us to write i?| = {2/mujQ){pQ — 
h{l - 0/(1 0) and i?2 = (2/mujl){pQ + h). 

In the superfluid phase all four elements of the com- 
pressibility matrix are equal. This structure results from 
the constraint p| = p|. In the fully polarized phase, 
the only nonzero element of the compressibility matrix 
is Thus in each region the compressibility matrix 

has rank 1, and there is only a single dynamical vari- 
able at each point in space. In the superfluid phase 
(r < we formalize this observation by taking the 

dot product (from the left) of the vector (1,-1) and 
Eq. The time derivative term vanishes and we are 
left with V • [(p|/m)V((5^| — Sp^)] = 0. A sufficient con- 
dition for this to be satisfied is 5p^{r) — Spi{r) = A, 
where A is independent of r (for r < Similarly, 
in the polarized phase (r > pi — and we get 

PjV • [{pi/m)V{Spi)] = 0. Again, a sufficient condition 
is 6pi{r > Ri) — s, where s is a constant. Motion of the 
boundary is captured by the time dynamics of s and A. 

To match the solutions in the two regions we note that 
at the edge of the minority cloud, the density is discon- 
tinuous. This implies that the compressibility matrix has 
a delta-function singularity. Including this singularity in 
Eq. Q, one sees that the quantity in the divergence, 
fa = VSpa/pa-, is discontinuous. By integrating Eq. 
in the neighborhood of one finds that this disconti- 
nuity is due solely to the discontinuity in pa-, and that 
both Spij and its first derivative are continuous. 

Equating the two expressions for Spi at the boundary, 
we have Sp^Ri) — s = X. Introducing w(r) = (5/X|(r) — 
5pi(Ri), we then have. 
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r<i?i, (6) 
r>Ri. (7) 



Subtracting Eq. © from 10 and setting r = Ri we get 
an ordinary differential equation for A, which is readily 
solved in terms of w{r = R^), allowing us to eliminate A 
from Eq. 10 and {T)) . Finding the compressibilities from 
the equation of states, we derive a closed set of equations 
for r > Ri and r < i?| respectively. 
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r^W^Sp^^iRl). Here R^ ^ R^^ - Rj, ^ - R^ 
and 5/if+/_(i?j^) is the fiuctuation of the majority species 
chemical potential on the normal/superfluid side of the 
boundary. Assuming Spi cx exp[iwt], and using a suitable 
change of variables, each of the equations in Eq. © 
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FIG. 1: Breathing mode frequencies of a spherical two com- 
ponent Fermi gas as a function of polarization P at unitar- 
ity. Solid line: From matching the analytic solutions at the 
boundary between superfluid and normal region in the two 
shell approximation. Dotted line: From numerical solutions 
of Eq. Q using the two-shell equation of state. Dashed line: 
Numerical solutions of Eq. in the BCS approximation 
(which predicts a three shell structure). 

can be converted into hypergeometric equations and the 
solutions are given in terms of Hypergeometric functions 
F{a, 6; c, u) as [l5| . 

6fx^ = CiF{a,p-r,,l-r^/R^^)r%„^{eA) (9) 
6m = -X + C2F{a,p-^y/6\'Yu„,{e,4>) (10) 

for r > Ri and r < i?| respectively. The parameters 
a = (l/2)[2 + I + V4 + 3e-hP + /], P = (l/2)[2 + I - 
V'i + 3e + P + 1], 7 = (2; + 3)/2, and 77 = a + /3 + 
1 — 7 are functions of magnetic quantum number I and 
the dimensionless mode frequencies e — oj'^/ojq. Here 
X = {2/3)[R^CiGi{R^)-r^C2G2{Ri)]R\nrn{0,(b). No- 
tice, we introduced two arbitrary constants Ci and C2. 
The functions Gi and G2 are defined as, Gi{Ri) = 
{a/h)F{a + 2,13 +2; 7^ + 2,1- Rl/R]) - {c/d)F{a +1,(3 + 
l;?7-M,l-i?2/i?2) and G2{Ri) = {a/ e)F{a + 2, (3+2;-^ + 
2, RyS^) + {c/{jS^)F{a + 1, /5 + 1; 7 + 1, Rf/S^) 

We match the two solutions and their derivatives at 
the boundary to get the frequencies lo for the breathing 
modes (/ = 0) at unitarity. The lowest energy breathing 
mode frequency is shown in Fig. ^ as a function of po- 
larization. We compare these results with those from the 
BCS equation of state. We also numerically solve Eq. ^ 
using the two-shell equation of state. The slight devia- 
tions of numerical and analytic results obtained using the 
two shell equation of state are numerical artifacts due to 
the smearing out of singularities in the compressibilities 
at the boundary. 

At both P = Q and P = 1 the whole cloud is in a 
single phase with an equation of state of the form n oc 
/i'^/^. Consequently, in both of these limits the cloud 



FIG. 2: Numerically calculated breathing mode frequencies 
in the BCS approximation. Dotted line: BCS regime where 
theory predicts a two/three shell structure; superfluid, par- 
tially polarized normal and fully polarized normal shells. 
Solid line: Crossover-BEC regime where theory predicts a 
two shell structure; superfluid and fully polarized normal 
shells. Dashed line: Deep BEC regime where theory predicts 
a two/three shell structure; superfluid, polarized superfluid, 
and fully polarized normal shells. These three graphs use the 
dimensionless parameter ^2mjj,o/h^as= -2.0, -1-2.0, -1-0.2 re- 
spectively. 

breaths with frequency to = 2ujo |^. For intermediate 
polarizations, the mismatch of the speed of sound causes 
a drop in the oscillation frequency. 

Away from unitarity we must rely on the mean-field 
calculations, as we know of no non-unitary analogy of 
Chevy's equation of state. In the BCS regime (os < 0), 
mean-field theory predicts a three-shell structure at low 
polarizations: an unpolarized superfluid core surrounded 
by partially polarized normal, and fully polarized nor- 
mal shells. At higher polarizations, the superfluid core 
is absent. Since in the absence of a superfluid re- 
gion the mean-field equation of state reduces to that of 
non-interacting particles, one finds that at these large 
polarizations the breathing mode frequency again be- 
comes u! — 2u!q (see FIG. EJ. In the deep BEC limit 
(flg > 0,kfas «C 1), the mean-field theory again pre- 
dicts a three-shell structure: an unpolarized superfluid 
core surrounded by a partially polarized superfluid shell 
and a fully polarized normal shell. At high polarizations 
the inner core will be absent. In this regime we do not 
see any dramatic signature of the disappearance of the 
central superfluid core. 

As seen in FIG. 12 the qualitative behavior of the 
polarization dependence on breathing modes has the 
same non-monotonic behavior in the entire BCS-BEC 
crossover region. 

In ref. we used a sum rule approach to find col- 
lective mode frequencies for unpolarized Fermi gases. 
Although the sum rule technique only provides upper 
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bounds on these frequencies, we found that the bound 
was very tight, and that metho d p roduced excellent 
agreement with both experiments |llL Il2l | and hydrody- 
namic theories |^ 0| . Repeating those calculations for 
the partially polarized gas, we find that the sum rules 
provide a much weaker upper bound in the present case. 
For example, using the two-shell equation of state at uni- 
tarity, the sum rule calculation finds no polarization de- 
pendance. It simply bounds lo < 2ll!q for all P. 

Although here we only consider spherically symmet- 
ric traps, we believe that the qualitative behavior of the 
polarization dependence of axial and radial breathing 
modes in an anisotropic traps will be similar to what 
we found. In particular the frequencies of these modes 
should be sensitive to the equation of state. We specu- 
late that in highly asymmetric traps, such as those used 
at Rice, surface tension effects may begin to play a role 
in the collective mode frequencies, though we have not 
calculated these effects. 
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